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We report improved measurements of the product branching fractions B(B + —► D°D* q"(2317)) x 
B(D*+( 2317) -> D+7t°) = (8.0i^±l.l±0.4)xl0“ 4 and B(B° -> D~D* s + (2317)) xB(L>*+ (2317) -A 
D+7r°) = (10.2++| ± 1.0±0.4) x 10 -4 , where the first errors are statistical, the second are systematic 
and the third are from D and Df branching fractions. In addition, we report negative results from 
a search for hypothesized neutral (z°) and doubly charged (z++) isospin partners of the D*g"(2317) 
and provide upper limits on the product branching fractions B(B° —> D° z°) x B(z° —> Dfn~) 
and B(B+ -A D~ z ++ ) x B{z++ ->• D+ 7 r+) that are more than an order of magnitude smaller than 
theoretical expectations for the hypotheses that the D* q" (2317) is a member of an isospin triplet. 
The analysis uses a 711 fb _1 data sample containing 772 million BR-meson pairs collected at the 
Y(4S) resonance in the Belle detector at the KEKB collider. 

PACS numbers: 12.39.Mk, 13.20.He, 14.40.Lb 


INTRODUCTION 

The (2317) meson, hereinafter referred to as the 
D*q , was first observed by BABAR as a narrow peak 
in the D + i r° invariant mass spectrum produced in in¬ 
clusive e+e~ —> D+ir°X annihilation processes 0,0, 


and confirmed by CLEO 0. Its production in the 73- 
meson decay processes B —> DD*^ was subsequently es¬ 
tablished by both Belle 0 and BABAR Q. (Here, B 
and D are used to denote B° and D~ or i3 + and D °.) 
Although it is generally considered to be the conventional 
I(J P ) = 0(0+) P -wave cs meson, its mass, M D ,+ = 
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2317.8±0.6 MeV [6, 7 ], is the same as the peak mass of its 
nonstrange counterpart, the 0 + P -wave cq (q = u or d ) 
Dq with mass Mjj* = 2318 ± 29 MeV [f|, in spite of the 
fact that the mass of the s quark is ~ 100 MeV above 
that of either of the q quarks. Potential-model Q and lat¬ 
tice QCD [ 9 ] calculations published prior to the BABAR 
discovery predicted that the 0 + P -wave cs meson mass 
would be well above the m^o + mx+ = 2358.6 MeV 
threshold and have a large partial decay width for the 
strong-interaction- allowed process D*f —> DK. The 
observation of a subthreshold mass has led to theoreti¬ 
cal speculation that the D* g~ is not a simple cs meson, 
but instead a DK molecule 0, a diquark-diantiquark 
state Em or some mixture of a cs core state with a DK 
molecule and/or a diquark-diantiquark [12]. 

A cs meson with mass below the 2358.6 MeV thresh¬ 
old would decay via the isospin-violating process D*f —>■ 
D + 7T° or the electromagnetic process D*f —► D * + 7 
and, thus, have a narrow natural width. This is con¬ 
sistent with experimental measurements, which have es¬ 
tablished a 95% C.L. upper limit on the total width of 
r D .+ < 3.8 MeV Q. The small width of the D is 
evidence for an I = 0 assignment. However, the CLEO 
experiment has established a stringent 90% C.L. upper 
limit on the partial width for D*f —» Df 7 decay Q: 


R(D*+) = 


r (d*+ Df7) 


F(D*+ D+ttO) 


< 0.059, 


(1) 


while studies that consider the D*f to be the cs chiral 


partner of the Df 13] predict values for R(D*f ) that are 
higher than the CLEO upper limit. Product branching 
fractions for B —> DD*f, D *q“ —> D+7r° have been 
measured by BABAR Sand Belle dj; the Particle Data 
Group (PDG) averages [6| of their results are: 


B(B+ D°D*+) x B(D*+ -»■ Df tt 0 ) = (7.3±?;f) x 10" 4 , 
B(B° -)■ D~D*+) x B(D*+ -j. D+tt 0 ) = (9.7± 4 ;°) x 10" 4 . 

Under the plausible assumption that B(D*f —> Df 7r°) ~ 

1, these measurements translate into the branching frac¬ 
tion ratios 


B(B+ D 0 D* s f) 
B{B+ -+ D°Df) 
B(B° -> D-D*+) 
B(B° D~Df) 


0 . 0811 q;o 2 i, 

0-i3±8;8g, 


which the authors of Refs. [14| and [15| note are well 
below expectations for a purely cs quark-antiquark state 
and an indication of some kind of multiquark content. 

The BABAR and Belle measurements for both B + and 
B° modes agree within errors, the biggest difference is 
1.5cr for the B° mode. In both cases, the measurements 
are based on event samples that are about 20% of the 
currently available data. Updated measurements based 


on the full data sets from both experiments would be 
useful. 

A report by Hayashigaki and Terasaki jl(| concluded 
that an I = 1 and I 3 = 0 assignment for the D*f can¬ 
not be ruled out and claimed, in fact, that an I = 1 
diquark-diantiquark interpretation is favored by some ex¬ 
isting data. If this were the case, doubly charged I 3 = 1 
(z ++ ) and neutral I 3 = —1 (z°) partners of the D*f 
with mass within ~ ±10 MeV of M n »+ should exist. 

^sO 

Since the z ++ and z° would be charmed mesons with 
/=1 and S = 1, they would necessarily have a minimal 
quark content of csud and csdu, respectively. Although 
a BABAR search for doubly charged and neutral part¬ 
ners of the D*f in inclusive e + e _ annihilation events sets 


95% C.L. upper limits on their production rates at 1.7% 
and 1.3%, respectively, of that for the D*f [17|, Terasaki 


has argued that these do not conclusively rule out their 
existence [lj|]. If the z ++ and z° mesons exist, isospin 
invariance ensures that the product branching fractions 
B(B -± D z++’°)_x B{ z++’° -± Dfn+’-) will be nearly 
equal to B(B -> DD* s f) x B{D* s f -> Dfn°). 

Here, we report measurements of B(B + —> D°D*f) x 
B(D*+ Df tt 0 ) and B(B° ^ D~D*+) x B(D*+ 
Dfn°) using a data sample that is more than 6 times 
larger than that used in previous results 0 and a 
search for doubly charged (z; ++ ) and neutral (z°) isospin 
partners of the D*f in the decay processes B + 


D~ z 


++ 


y + + 


Df 7r + and B° 


Df 7 


The results are based on the full Belle T(4S) data sam¬ 
ple (711 fb _1 ) that contains 772 million HH-meson pairs 
produced at a center-of-mass system (cms) energy of 
y/s = 10.58 GeV and collected in the Belle detector at 
the KEKB energy-asymmetric e + e~ collider jl9| . 


DETECTOR DESCRIPTION 

The Belle detector is a large-solid-angle magnetic spec¬ 
trometer that consists of a silicon vertex detector, a 
50-layer cylindrical drift chamber, an array of aerogel 
threshold Cherenkov counters, a barrel-like arrangement 
of time-of-flight scintillation counters, and an electromag¬ 
netic calorimeter comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 
1.5 T magnetic field. An iron flux-return located out¬ 
side of the coil is instrumented to detect mesons and 
to identify muons. The detector is described in detail 
elsewhere po} . 


EVENT SELECTION 

We reconstruct Df mesons via their tt + K + K~ de¬ 
cay mode, which has a branching fraction of B D + = 
(5.39 ± 0.21)%, D~ mesons via the A' + 7r”7r _ decay 
mode [B d - = (9.13 ± 0.19)%] and 5° mesons via the 
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K+n~ [Bkk = (3.88±0.05)%1 and K+n+TT~n- [lB K37r = 
(8.08 ± 0.20)%] decay modes fj]. 

For all charged particles, we require dr < 0.7 cm and 
\dz\ < 3.0 cm, where dr and dz are the track’s dis¬ 
tances of closest approach to the run-dependent mean 
interaction point transverse to and parallel to the e + 
beam direction, respectively. Charged-particle identifi¬ 
cation is accomplished by combining information from 
different detector subsystems to form likelihood ratios, 
Lr/v = Lk/(Lk + I/7I-), where Lk {L n ) is the likelihood 
of the kaon (pion) [2l[. A charged track is classified as 
a kaon (pion) if L K ^ t u/k) > 0-5, with both the muon 
likelihood ratio and electron likelihood smaller than 0.95. 
For B° —> D~D*q decay, the kaon and pion identification 
efficiencies both exceed 95%. We reconstruct ^ r° mesons 
via their 7r° —> 77 decay mode using 7 candidates with 
E-y > 30 MeV and 77 combinations that satisfy a one- 
constraint (1C) kinematic fit to m n 0 with y 2 < 6.0. In 
addition, we require |M 77 — m n o| < 12 MeV and the 7r° 
three-momentum in the e + e _ ems p™ s < 1.9 GeV. 

Candidate D mesons are required to have a Knit 
(n = 1 to 3) invariant mass in the range \Mk U 7T — mo\ < 
2.5 a of the observed peak mass, where a is the width 
from a Gaussian fit to the Krm invariant mass peak; 
D+ candidates are required to be in the mass interval 
I M k+ K-K+ — m D+\ < 2.5ct. Here, the values of cr range 
from 4.6 MeV to 5.5 MeV. 

Candidate B -» DD*q decays are identified by i) the 
ems energy difference A E = E™ s — ii) the beam- 

energy constrained mass M bc = vA-^Wm) 2 — () 2 > 
and in) the Dfn 0 invariant mass. Here -E£“ s m is the 
ems beam energy and Fig 118 and are the total 

ems energy and three-momentum of the particles form¬ 
ing the DD*q combination. We select events with 
Mb c > 5.20 GeV, -0.12 GeV < A E < 0.1 GeV and 
2.228 GeV < M D +^ 0 < 2.418 GeV for three-dimensional 
fitting, and define signal regions as |M bc — tub\ < 
0.007 GeV, -0.033 GeV < A E < 0.030 GeV and 
\M D + n0 — 2. 3178 GeV| < 0.015 GeV. For candidate B -> 
Dz ++ (z°) decays, the 7r° is replaced by a 7r + (n~) and 
the A E signal region is compressed to | AE\ < 0.023 GeV. 
These intervals correspond approximately to ±2.5er win¬ 
dows around the central values for each variable. 

To reduce background from e + e _ —> qq continuum 
processes, where q = u,d,s,c, we require the follow¬ 
ing: R 2 < 0.3, where R 2 is the normalized second Fox- 
Wolfram moment [22j]; |cos(?b| < 0.8, where 9b is the 
polar angle of the candidate S-meson direction in the 
ems; and | cos0thrs| < 0.8, where 0thrB is the ems angle 
between the thrust axis of the B candidate and that of 
the remaining unused tracks in the event. These require¬ 
ments reject 14% of B° —$■ D~D*jf~ signal and 45% of qq 
continuum. 


MC SIMULATION 

We use Monte Carlo (MC) simulation to optimize se¬ 
lection criteria, determine acceptance and study multiple 
candidates per event [23|. We generate signal MC for 
each process under investigation using PDG values 6;] 
for subdecay branching fractions and setting B{D*q —> 
U+7T°) and B( z ++, ° —> U+7r +,_ ) = 1. In addition, we 
use a generic BB MC sample with about 3 times the inte¬ 
grated luminosity of the actual data sample to investigate 
possible peaking backgrounds. The simulated events are 
processed through the same reconstruction and selection 
codes that are used for the real data. 


MULTIPLE CANDIDATES 

The D*jf~ —► .D+7T 0 mode is plagued by a large frac¬ 
tion of events with multiple candidates. The numbers of 
events with multiple entries in the full fitted region are 
summarized in Table |I] Since the MC samples reproduce 
the data reasonably well, we use the MC as a guide for 
methods to reduce the multiple candidates. 


TABLE I: Fractions of multiple candidate events in data and 
MC. 


Sample 

->• D-ET+ 

D — > K 7T7T 

B+ ->• D U D* S + 
D° -»• Kn 

B+ ->• D U D*+ 
D° -> K3n 

Sig. MC 

70% 

45% 

70% 

BB MC 

69% 

39% 

69% 

Data 

68% 

39% 

69% 


For the D~ and D° —> 

modes, about two thirds of the multiple candidates are 
low-energy photons forming multiple 7r° —> 77 combina¬ 
tions and one third are multiple charged pions in the D 
candidate. For the D° —> K + n~ mode, essentially all of 
the multiple candidates are associated with the 7r° —> 77 
reconstruction. 

We use the 77 energy asymmetry, E asym = (E\ — 
E 2 )/(E 1 + E 2 ), where Ei (E 2 ) is the higher-(lower- 
energy photon of the 77 pair, to select 7r° candidates. 
The left panel of Fig. [I] shows the -E a sym distribution for 
correctly assigned 77 pairs in signal MC events; the right 
panel in the same figure shows the same distribution for 
incorrectly assigned combinations. Here, the events are 
required to be in the M bc and A E signal regions. Accord¬ 
ing to MC studies, the strong peak near E asym ~ 0.85 
in the incorrect-assignment plot is mostly due to beam- 
produced background photons. Figure [2] shows the cor¬ 
responding x 2 distributions from the 7r° —> 77 kinematic 
fits. To reduce the 7-associated multiple candidates while 
minimizing loss of signal efficiency, we require that pho¬ 
tons in the energy interval 30 MeV < 2£ 7 < 40 MeV have 
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X 2 < 0.5 for the 1C fit or E asym < 0.7. For remaining 
events with multiple 7 candidates, we select the combi¬ 
nation with the smallest E asym value. For multiple D 
(Df ) candidates, we select the track combination with 
invariant mass closest to the PDG value for mu (m D +). 



FIG. 1: The E asym distributions for signal MC events for 
correctly (left) and incorrectly (right) assigned photons. 



peaking background added and the solid blue curve in¬ 
cludes the signal. 




FIG. 2: The y 2 distributions from the n° —> 77 fit for 
signal MC events for correctly (left) and incorrectly (right) 
assigned photons. 


DDg+ EFFICIENCIES 


We determine event yields from unbinned three- 
dimensional likelihood fits [M bc vs M(Df n 0 ) vs A E] 
to the selected data using a bifurcated Gaussian function 
for the Mb c signal probability density function (PDF) 
and an ARGUS function |24| multiplied by a second- 
order Chebyshev polynomial for the M bc combinatorial- 
background PDF. For A E, we use a Crystal Ball func¬ 


tion [25] for the signal PDF and a third-order Cheby¬ 
shev polynomial for the combinatorial-background PDF. 
For M(Dfir°), we use a Gaussian function for the sig¬ 
nal PDF and a third-order Chebyshev polynomial for the 
combinatorial-background PDF. 

In the generic BB MC samples, there is background 
that peaks in M bc and A E [but not M(Dfn 0 )] mostly 
coming from three-body B —y Dtt q D£ decays. This 
background is modeled in the fits by M bc and A E signal 
functions and a linear function for 

As an example, we show fit results for the B° —> 
D~D* s q“ signal MC sample in the upper part of Fig. [3] 
The lower part of Fig. [3] shows the results from fits to 
the generic MC sample. In these figures and subsequent 
plots in this report, the red short-dashed curve is the 
fitted background; the green long-dashed curve has the 


FIG. 3: Top: The M bc (left), M(L>+tt 0 ) (center) and A E 
(right) distributions for the B° —y D~D*q signal MC events 
with the results of the fit superimposed. The events in each 
distribution are in the signal regions of the two quantities not 
being plotted. Bottom: The corresponding distributions for 
the generic MC event sample (~ 3 times the data). (See text 
for curves.) 

The detection efficiencies determined from the signal 
MC events that survive the application of the multiple 
event selection requirements are listed in Table El 


TABLE II: The MC-determined B —y DD*q efficiencies. 



D —3 K + 7T 7T 

B + D U D*+ 
D° ->• A'+tU 

B+ D U D*J 

D° —y I\ + 7T + 7T 7r 

Agen 

266230 

266230 

266230 

N &t 

7022 ± 90 

8575 ± 97 

4839 ± 72 

effic. 

(2.64 ± 0.03)% 

(3.22 ± 0.04)% 

(1.82 ± 0.03)% 


B -y DDg ^; D* s + -y .D^tt 0 RESULTS 

1 ) b° d-d*+, d:+ d+ tt° 

We determine the number of B° —y D~ D*^; D*q —> 
Dfn° signal events in the data by applying the three- 
dimensional fit described above to the selected D = D~ 
candidates. In this fit, the rms widths of the Mb c , 
M(Df 7 T°) and A E signal functions are kept fixed at their 
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MC-determined values. Figure [4] shows the results of 
the fit, which returns a signal yield of N evt = 102.6^2 o 
events. The fitted peaking background yield is consis¬ 
tent with zero: 7.7± 13.6 events. The signal significance, 
determined as the square root of twice the difference of 
log-likelihood values from fits with and without a signal 
term, is 9.9 a. 



FIG. 4: The Mbc (left), M(Dfn°) (center) and A E (right) 
distributions for projections of the B° —¥ D~ D*f candidate 
events that are in the signal regions of the two quantities not 
being plotted. The results of the fit described in the text are 
superimposed. (See text for curves.) 

We determine the product branching fraction from the 
relation 

B(B° D~D*+) x B(D*+ Dfn°) (2) 

-^evt 

^BBVd-D+^ d -^D+ 

where N BB = (772±11) x 10 6 is the number of BB events 
in the data sample and V D - D + is the MC-determined 
detection efficiency for this channel (see Table m • The 
result is 

B(B° D~D*+) x B(D*+ Dfn°) (3) 

= (10.21^ ±1.0 ±0.4) x 10~ 4 , 

where (and elsewhere in this report) the first error is sta¬ 
tistical, the second is the systematic error (discussed be¬ 
low), and the third reflects the errors on the PDG branch¬ 
ing fractions of the D~ and Df mesons [if. This result 
agrees well with the average of the BABAR and previous 
Belle measurements mentioned above with a substantial 
improvement in precision. 

2 ) B+ -»■ D 0 D;+, D*+ D+ir° 

The top plots of Fig. [5] show the Afbc, M(D+7r°) and 
A E distributions of the B + —> D°D*f, D*q“ — > Dfir 0 , 
D° —> K + tt~ candidates. Here, in addition to the rms 
widths, we fix the Mb c and A E peak positions. The fit 
results are 38.9+g2 signal events and 12.6^ 22 7 6 peaking 
background events. An application of the equivalent of 


Eq. (]2j) to this mode results in the product branching 
fraction 

B{B+ D°D*+) x B(D*+ -+ Df tt 0 ) (4) 
= (7.5lf g ± 0.7 ± 0.3) x 10 -4 , 

which is in good agreement with the PDG average of 
previous measurements but with a smaller error. 

The bottom plots of Fig. [5] show the Mb c , M(Dfn°) 
and A E distributions of the B + — > D°D* g“, D*f —> 
D+tt°, 5° -> K+TT+ 7 r tt candidates. Here again, in 
addition to the rms widths, we fix the M^ c and A E peak 
positions. The fit results are 52.4+ 42 'g signal events and 
99.0ig;§ peaking background events. An application of 
the equivalent of Eq. © to this mode results in the prod¬ 
uct branching fraction 

B(B+ -j. D°D*+) x B(D*+ -j. Dfn°) (5) 
= (8.6l 2 ;g ± 1.1 ±0.4) x 10“ 4 , 

which is in good agreement with the result for the D° —> 
K + tt~ mode and the PDG average of previous measure¬ 
ments and with a comparable error. 

The weighted average of the two measurements is 

B(B+ -j. D°D*+) x B(D*+ Dfn°) (6) 
= (8.0±i;l ± 1.1 ±0.4) x 10" 4 , 

where near-complete correlation of the systematic errors 
for the two measurements is taken into account. 

As a consistency check, we apply a simultaneous fit 
to the two modes, where we find a total signal yield 
of 91.9lf^ with a statistical significance of 5.9o\ The 
peaking background yield is 148. 5 + 2^5 even t s - The sig¬ 
nal yield from the simultaneous fit is consistent with the 
sum of individual fits, while the number of peaking back¬ 
ground events is marginally higher. The product branch¬ 
ing fraction obtained using the simultaneous fit is 

B(B+ -j. D°D*+) x B(D*+ -j. Dfn°) (7) 
= (8.11^3 ± 1.1 ± 0.3) x 10" 4 , 

in good agreement with the result from the weighted av¬ 
erage of results for each mode. 

3) Systematic errors 

Systematic errors include the errors on N BB and the 
D and Df secondary branching fractions, MC statis¬ 
tics and model dependence, MC-data differences in par¬ 
ticle identification, charged-particle tracking, 7r° identi¬ 
fication, and the choice of the fitting model. The error 
on N bb is 1.4% and the secondary branching fraction 
relative errors are the PDG values: D + —> K~7 r + 7r + 
(2.0%); D° K~tt+ (1.3%); D° R-tt+tt+tt- (2.6%); 

















7 




5.2 5.22 5.24 5.26 5.28 5.3 2.25 2.3 2.35 2.4 -0.1 -0.05 0 0.05 0.1 


M bc (GeV) M(D* jt°) (GeV) AE(GeV) 

FIG. 5: Top: The M bc (left), M(D+ n°) (center) and A E 
(right) distributions for the B + —> D° D*f candidate events 
for the D° —> I\ + n~ subdecay mode, with the results of the 
fit superimposed. The events in each distribution are in the 
signal regions of the two quantities not being plotted. Bot¬ 
tom: The corresponding distributions for D° —» K + tv + tt~'k~ 
decays. (See text for curves.) 

Df —> K + K~n + (3.9%). The MC model dependence 
is evaluated by varying the Df —> cf>ir + component of 
Df —> K + K~tt + decays between extreme limits and 
changing the phase-space distributions for the multibody 
D-meson decay modes. We use various control samples 
to determine MC-data efficiency differences that are com¬ 
mon to many Belle analyses to evaluate systematic er¬ 
rors associated with kaon (pion) identification of 1.1% 
per track (1.2% per track), charged particle tracking of 
0.35% per track, and i r° detection of 4.0%. 

The dependence on the fitting model is estimated from 
changes observed by redoing the fits with each parame¬ 
ter fixed at ±1 ct from its best-fit value. The systematic 
errors from each source, listed in Table cni are summed 
in quadrature to get the final value. 


SEARCH FOR z++ -D+7T+ AND z° -j- D+n~ 

We look for z ++ — > Df tt + and z° —»• Df tt~ signals in 
the B + —> D~D+tt + and B° —> D°Df tt~ decay chan¬ 
nels by applying the selection criteria discussed above 
with the replacement of the selected 7 r° with a 7r + (for 
z ++ ) or 7T _ (for z°). Here, for events with multiple 
D and/or Df track combinations, we select those with 
a measured invariant mass closest to the corresponding 


TABLE III: Summary of relative systematic error sources (in 
percent). 



B u -> D-Dlf 

B+ —> D u D*f 

B+ D u D*f 


D —> Knn 

D° -S- Ktv 

D° K3n 

D&D+ BFs 

4.4 

4.1 

4.7 

Nbb 

1.4 

1.4 

1.4 

MC model dep. 

3.6 

2.3 

5.9 

MC stat. 

1.2 

1.0 

1.4 

Particle ID 

6.9 

5.2 

8.4 

Tracking 

2.1 

1.8 

2.5 

Fit params. 

4.4 

5.8 

4.7 

n° 

4.0 

4.0 

4.0 

Quad, sum 

10.2 

9.4 

12.4 


PDG values. For z ++ signal MC, the number of remain¬ 
ing events with multiple candidates is 11.2% over the 
full three-dimensional range of the likelihood fit; for z°, 
fewer than 0.1% of the remaining events have multiple 
candidates. 

1) Peaking backgrounds from generic MC samples 

We check for possible peaking backgrounds leaking into 
the signal using a sample of simulated generic B-meson 
decay events (with no z ++ or z° signals) with a luminos¬ 
ity that corresponds to 3 times the number of B decays 
in the data. The top plots of Fig. [6] show the results of 
applying the three-dimensional fit to selected D~ Df n + 
MC events. Here, the signal yield is zero with a posi¬ 
tive error of 7.1 events. The peaking background yield is 
544±41 events. The middle (bottom) plots of Fig.[G]show 
the results of the three-dimensional fits to the generic MC 
for the D° —> K + n~ ( D° —> I\ + tt + t:~t:~) channel in the 
selected B —> z°D° samples. No background processes 
are found that produce a spurious signal; the signal yields 
are also zero for both D° modes with positive errors of 
2.1 and 9.9 events for the K + tt~ and Jv + 7r + 7r~7r - modes, 
respectively. The M^, C -AE peaking background yields for 
these modes are 169 ±22 and 229^32 events, respectively. 

2) Mass-dependent efficiency 

Since the z ++ and z° are hypothesized to be isospin 
partners of the B*q", their masses are expected to lie 
somewhere within a ±10 MeV mass region of m D *+ = 
2317.8 ± 0.6 MeV. In order to be certain that we cover 
all reasonably plausible mass values, we scan for z ++ and 
z° signals in 13 adjacent mass bins, each 5 MeV wide, 
covering a ±32.5 MeV interval centered on 2317.8 MeV. 

To account for a possible mass dependence of the detec¬ 
tion efficiency, we generate z ++ and z° signal MC events 






























































are summarized in Table m As examples, we show 
the fit results for the mass bin centered at M(Df ir) = 
2317.8 MeV for the z ++ (z°) search in the top (bottom) 
plots of Fig. [3 None of the fits returns a positive z ++ 
or z° signal with a statistical significance of more than 
1.3cr. The determination of the Bayesian 90% credibility 
level upper limits |26| on the event yields and product 
branching fractions is described below. 


TABLE IV: Product branching fraction upper limits £5+ for 
B(B+ ( B° ) -)• D~ ( D°) Zi ) X B{Zi -> D+tt ) ( Zl = z++ and 
Z 2 = z°), for Zi masses between 2285.3 MeV and 2350.3 MeV. 
Here AM = M c t r — m n *+, where M c t r is the center of the 

sO 

5 MeV mass window allowed for the fit, and Ni is the upper 
limit including systematic errors. 


AM 


/yUG 

R ul 

D ++ 

A 0 nt 

V 0 UL 

So" 

MeV 



(to- 4 ) 



(io- 4 ) 

-30 

4.o±S:S 

16.3 

0.52 

-13.7 ±6.2 

10.5 

0.34 

-25 

4.ii 5 s ;^ 

16.3 

0.52 

c c+7.y 

21.2 

0.69 

-20 

-s.3±t;l 

9.8 

0.32 

5.8±|;° 

21.5 

0.69 

-15 

-10.31^ 

8.0 

0.25 

2.7 ±8.3 

20.1 

0.65 

-10 

-10.2 ±3.5 

7.9 

0.25 

^• U -8.4 

20.4 

0.66 

-5 

-8.8 ±3.2 

8.1 

0.25 

4.1 ± 7.4 

20.4 

0.66 

0 

-9.3 ±3.0 

8.4 

0.27 

o 1+7.8 

19.8 

0.64 

5 

-9 3+ 4 - 5 

8.5 

0.28 

1 y+10.3 
-*-• ‘ -6.1 

16.0 

0.52 

10 

4 6 +5 - 6 

16.2 

0.51 

-5.41S1 

13.4 

0.44 

15 

6.4 ±5.0 

17.8 

0.57 

— 5-4+J 

13.3 

0.43 

20 

6.oi 5 s .? 

17.6 

0.56 

q q+H.5 

14.3 

0.47 

25 

3-0±«;g 

15.8 

0.50 

5-7±5;g 

20.6 

0.67 

30 

q 4+5-7 

15.8 

0.50 

5.611;? 

20.0 

0.65 


FIG. 6: The Mb c (left), M(D+7r) (center) and A E (right) 
distributions for generic-MC events that pass the D~Dfn + 

(top), D°Df 7r~, D° -*■ K + n~ (middle) and D°Df n~, D° -)• 

A' + 7r + 7r - 7r^ (bottom) channels. The curves are the results of 
fits described in the text. 

4) Systematic errors for z ++ and z° searches 


with z masses in the full range of the scan. The effi¬ 
ciencies, determined from fits to the selected events from 
each MC sample, are independent of mass to within the 
~2.5% MC statistical errors. For the z ++ search, the 
average efficiency is (8.3 ± 0.1)%. For the z° search, the 
average efficiency is (9.2 ± 0.1)% for the D° —> K + ir~ 
mode and (4.1 ± 0.1)% for D° — > K + n + n~n~. 

3) Fits to the AF(D+7r + ’ — ) spectra 

We apply a sequence of 13 three-dimensional fits to the 
data using a Gaussian signal function with width fixed at 
the MC-determined Df 7+ mass resolution (cr=4.6 MeV) 
to represent the z ++ (z°) with a peak mass restricted 
to 5 MeV-wide windows covering a total mass range of 
±32.5 MeV about m n .+ = 2317.8 MeV. The results of 
these fits for the z ++ —»• Df 7r + and z° —> Df tt~ searches 


Systematic errors are evaluated using the same meth¬ 
ods that are used for the branching fraction mea¬ 
surement described above, with the 7r°-associated error 
replaced by the error on the additional charged pion. 
For this, the nominal 0.35% tracking error is assigned to 
p > 200 MeV tracks. However, 5% of the relevant pions 
for the z° have p < 200 MeV with an associated error 
of 5%. Here, a weighted average is used and the total 
tracking uncertainty increases to 3.8%. For the system¬ 
atic error associated with multiple candidates, we per¬ 
form a multiple-candidate-free z ++ scan where we use 
the smallest A E to select the best candidate and a two- 
dimensional fit [Mbc and M(Dfn + )] to measure signal 
yields. From the differences between the results of the 
two methods, we determine a systematic error from this 
source of 2.2%. For other sources of error, we use the 
results listed in Table Hill The resulting errors are 11.4% 
for the z ++ search and 16.6% for the z° search. 
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FIG. 7: The M^ c (left), M(D^n) (center) and A E (right) 
distributions for selected B + —y D~D+n + (top) and B° —y 
D°Dfn~ (bottom) event candidates for the fit with the sig¬ 
nal peak mass restricted to a 5 MeV region centered at 
M(Df 7r) = 2317.8 MeV. In the lower plots, the D° —y K + n~ 
and A' + 7r + 7r _ 7r _ samples are combined. (See text for curves.) 


5) Upper limit determination 


where the notation follows that of Eq. m and rj ++ is 
the MC-determined efficiency. For the z° search, where 
there is no evidence for the signal either, we use the same 
relation with B D ~r] ++ replaced by BkttVktt + Bk 3 tvVk 3 -k, 
where t]k-k (vK 3 n) is the efficiency for the D° —y K + tt~ 
mode. The resulting product branching 
fraction upper limits, listed in Table HVl are all more than 
an order of magnitude lower than the measured values for 
the DD*q final states. This is in strong contradiction to 
expectations for the hypothesis that the D*q is a member 
of an isospin triplet [16]. 


SUMMARY 

We reported measurements of the product branch¬ 
ing fractions B{B+ —y D°D*+ ) x B{D*+ —y D+tt°) = 
(8.0±1;1 ± 1.1 ± 0.4) x 10" 4 and B(B° -> D~D*+) x 
B{D*+ -y D+tt°) = (10.2±i;| ± 1.0 ± 0.4) x 10" 4 . Here, 
the first errors are statistical, the second are systematic 
and the third are from D and Df branching fractions. 
These values agree with the existing PDG world average 
values Q, significantly improve upon their precision, and 
supersede those of Ref. 0 • In addition, we reported nega¬ 
tive results on a search for hypothesized doubly charged 
and neutral isospin partners of the and provided 

upper limits on the product branching fractions that are 
more than an order of magnitude smaller than the theo¬ 
retical predictions of Hayashigaki and Terasaki 0. 


We use a Bayesian method to convert the fitted results 
to upper limits on the total number of signal events. To 
account for the systematic uncertainties, the likelihood 
distributions from the z ++ (z°), fits are convolved with 
a Gaussian with (j syst = 0.114 (0.166) x j., where N^ t 
is determined from 

;v UL 

/ - C(n sig )dn sig = Q g (g) 

Jo £(n s i g )dn s i g 


The Gaussian width is cr sys t = 1.1 (3.1) events for the 
2317.8 MeV mass bin of the z ++ (z°) scan; the widths 
for the other mass bins are similar. The corresponding 
upper limits, iV UL , are determined from the relation 

jV UL 

Jo B(n s j g ) <^) Q(n s i g )dn s i g _ q g 
/ 0 +o ° £(n sig ) (g) G(n s i g )dn s i g 

and in all cases differ from N^ t by less than one event. 
The resulting values of IV UL are listed in Table m 
For the z ++ search, we determine upper limits on the 
product branching fractions U++ = B(B —y D~ z ++ ) x 
B{ z ++ —y Df 7r + ) from the relation 


UL _ 
++ ~ 


^BB^Dt^ D ~d++ 


( 10 ) 
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